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ABSTRACT 



We report the detection of a double planetary system around HD 140718 as well as the discovery of two long period and massive 
planets orbiting HD 171238 and HD 204313. Those discoveries were made with the CORALIE Echelle spectrograph mounted on the 
1.2-m Euler Swiss telescope located at La Silla Observatory, Chile. The planetary system orbiting the nearby G9 dwarf HD 147018 is 
composed of an eccentric inner planet (e=0.47) with twice the mass of Jupiter (2.1 Mj up ) and with an orbital period of 44.24 days. The 
outer planet is even more massive (6.6 Mj up ) with a slightly eccentric orbit (e=0.13) and a period of 1008 days. The planet orbiting 
HD 171238 has a minimum mass of 2.6 Mj up , a period of 1523 days and an eccentricity of 0.40. It orbits a G8 dwarfs at 2.5 AU. The 
last planet, HD 204313 b, is a 4.0 Mj up -planet with a period of 5.3 years and has a low eccentricity (e = 0.13). It orbits a G5 dwarfs at 
3. 1 AU. The three parent stars are metal rich, which further strengthened the case that massive planets tend to form around metal rich 
stars. 
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1. Introduction 

The CORALIE radial velocity planet-search program has been 
ongoing for more than 10 years (start date: June 1998) at the 1.2- 

Q\ meter Swiss telescope located at La Silla Observatory, Chile. It is 

r**"* a volume limited planet search survey that contains all Hipparcos 
main sequence stars from F8 down to K0 within 50 pc and has a 
color-dependant distance limit for later type stars down to M0 

QO (Udry et al. 2000). The fainter targets do not exceed V=10. 

f^*> Among the 1647 stars surveyed with CORALIE, 40 percent of 
them are measured with a radial velocity accuracy of 6 ms _1 or 
better and 90 percent of the sample is monitored with an accu- 
^ racy better than 10 ms _I . The remaining 10 percent of the sam- 

• <-h pie have measured at a lower accuracy due to the lower signal 
to noise ratio and/or to the the widening of the cross-correlation 
5-j function induced by large stellar rotation velocities. CORALIE 
went through a major hardware upgrade in June 2007 to in- 
crease the overall efficiency whith an net gain in magnitude of 
1.5 which allow to survey the fainter part of the sample with an 
accuracy of 5-6 ms _1 . 

So far, CORALIE has allowed the detection (or has con- 
tributed to the detection) of 54 extra-solar planet candidates. 



Send offprint requests to: Damien Segransan, e-mail: 
Damien. Segransan@unige . ch 

* Based on observations collected at the ESO La Silla Observatory 
with the CORALIE Echelle spectrograph mounted on the Swiss tele- 
scope. 
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This substantial contribution together with discoveries from var- 
ious other programmes have provided a sample - as of today - 
of close to 350 exoplanets that now permit to point out interest- 
ing statistical constraints for the planet formation and evolution 
scenarios (see e.g. Ida & Lin 2004; Marcy et al. 2005; Udry 
et al. 2007; Udry & Santos 2007; Mordasini et al. 2008, and 
references therein for reviews on different aspects of the orbital- 
element distributions or primary star properties). 

Concerning the upper part of the planetary mass distribution, 
it should be noted that a wealth of massive planets have been and 
is beeing discovered. For instance, 142 of the 350 discovered 
planets are more massive than 1 Mj up with periods larger than 
50 days. Their eccentricity distribution is similar to the binaries 
and only 18 of them are in multiple giant planet systems. Most 
of their parent stars are metal rich pointing toward the existence 
of disk enriched in heavy elements. Interesting enough, the mas- 
sive planet around the metal deficient star HD 33636 turns out 
to be a very low mass stars as shown by FGS/HST astrometric 
observations presented in (Bean et al. 2007). Although not sta- 
tistically significant, this result show how important astrometric 
observations are to clarify the nature of those massive planets. 

In this paper we report the discovery of a planetary system 
composed of two massive planets around HD 147018 and of 
two additional long period and massive planets in orbit around 
HD 171238, and HD 204313. Three of those planets come in 
addition to the 54 known planets with periods longer than 
1000 days and masses larger than 1 Mj up . The paper is orga- 
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nized as follows. In the second section, we discuss the host stars 
properties. The third section describes the instrumental upgrade 
of CORALIE, the resulting radial-velocity measurements and the 
orbital solutions. In section 4, we provide some concluding re- 
marks. 



2. Stellar characteristics 

Effective temperatures, gravities and metallicities are derived us- 
ing the spectroscopic analysis of Santos et al. (2000) while the 
v. sin (z) is computed using Santos et al. (2002)'s calibration of 
CORALIE's Cross-Correlation Function (CCF). We also used the 
improved Hipparcos astrometric parallaxes re-derived by van 
Leeuwen (2007) to determine the the V-band magnitude using 
the apparent visual magnitude from Hipparcos (ESA 1997). 
Metallicities, together with the effective temperatures and abso- 
lute V-band magnitudes are used to estimate basic stellar param- 
eters (ages, masses, radii and log g) using theoretical isochrones 
from Girardi et al. (2000) and a Bayesian estimation method 
described in da Silva et al. (2006). The web interface for the 
Bayesian estimation of stellar parameters, called PARAM 1.0 
can be found at http://stev.oapd.inaf.it/cgi-bin/param. Resulting 
Stellar parameters are listed in Table 1 . 

2.1. HD 14701 8 (HIP 80250) 

HD 147018 is a G9 dwarf with an astrometric parallax of n — 
23.28 ±0.86 mas and an apparent V band magnitude of V= 8.30. 
Our spectral analysis results in an effective temperature of T e ff = 
5441 ± 55 K and a stellar metallicity of [Fe/H] = 0.10 ± 0.07. 
Using theoretical isochrones, we finally derived a mass of = 
0.927 ± 0.031 M with an age of 6.36+4.33 Gyr. 

2.2. HD 1 71238 (HIP 91085) 

HD 171238 is a G8 dwarf with an astrometric parallax of n — 
19.89 + 1.15 mas and an apparent V band magnitude of V= 
8.61. Our spectral analysis results in an effective temperature 
of T eff = 5467 ± 55 K and a stellar metallicity of [Fe/H] = 
0.17 ± 0.07. Using theoretical isochrones, we derive a mass of 

= 0.943 ± 0.033 M Q with an age of 4.92+4.1 1 Gyr. 
Eventhough Hipparcos photometry is relatively stable, with a 
scatter of 0.018 mag in the visible, it should be noted that the 
star is listed in the General Catalogue of Variable Stars (Samus 
et al. 2009) as a BY Draconis-type variable. Such variable stars 
present photometric variability - induced by spot coverage or 
by chromospheric activity - up to 0.5 magnitude in the visible 
on time scales ranging from a fraction of a day to 120 days 
that are likely to affect the velocities. We were not able to de- 
rive a value of the log f^^) index since the star is too faint to 
conduct a proper spectral analysis with CORALIE, but we were 
able to averaged 35 CORALIE spectra taken with the Thorium 
lamp. As shown on Fig.l, a clear Ca II re-emission is seen at 
A = 3933.66 A revealing the presence of a significant chromo- 
spheric activity possibly induced by stellar spots or plagues. 

2.3. HD 20431 3 (HIP 106006) 

HD 2043 13 is a G5 dwarf with an astrometric parallax of n — 
21.11 + 0.62 mas and an apparent V band magnitude of V= 7.99. 
The spectral analysis results in an effective temperature of T e ff = 
5767 ± 17 K and a stellar metallicity of [Fe/H] = 0.18 + 0.02. 




3931 3932 3933 3934 3935 3936 3937 

Wavelength [Angstrom] 

Fig.l. Ca II H emission region for HD 171238. The large re- 
emission at the bottom of the Ca II H absoption line at A = 
3933.66 A is an indicator of chromospheric activity. This region 
is contaminated by a few thin thorium emission lines. Thorium 
pollution is too important to display the Ca II K region (around 
A = 3968.47 A). 

Using theoretical isochrones, we derive a mass of M+ = 1 .045 ± 
0.033 M Q with an age of 3.38+2.58 Gyr. 

Table 1. Observed and inferred stellar parameters of the plan- 
ets host's stars presented in this paper. (1) : Parameter derived 
from Girardi et al. (2000) models. (2) : Parameter derived using 
CORALIE CCF. 



Parameters HD 147018 HD 171238 HD204313 



3. Radial velocities and orbital solutions 

3.1. CORALIE upgrade 

Triggered by the interest to carry out spectroscopic follow-up 
on transit candidates fainter than the stars surveyed in our main 
planet search programme (V<10), we decided to to improve the 
overall efficiency of the instrument by upgrading CORALIE in 
June 2007. The fibre link and the cross-disperser optics have 
been removed and replaced by a new design. The double scram- 
bler has also been removed and the grism/prism cross-disperser 
component replaced by a series of 4 Schott F2 prisms of 32 deg 
angle each. The net outcome of this new design is to maintain 
the spectral range from 381 to 681 nm but with a large effi- 
ciency gain of about a factor of 6 (8 below 420 nm) and a spec- 
troscopic resolution of 55 000-60 000 (increased by 10-20%). 
Those hardware modifications have, however, affected the in- 
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strumental zero point with radial velocity offsets that could reach 
up to 20 ms _1 , depending on the target spectral type. For this 
reason, we decides to refer to the orginal CORALIE as CORALIE- 
98 and to the upgraded one as CORALIE-07. The overall instru- 
mental precision was not affected by the upgrate and stays at the 
5 ms~' level. 

The direct consequence - on our main planet search sur- 
vey - of the upgrade to CORALIE-07 is to to increase the 
efficiency of the instrument on bright targets (V<8.5) and to 
improve the radial velocity accuracy on the fainter part of our 
sample (28% of the sample, V=8.5-10). Those stars are now 
monitored with a long term radial velocity accuracy of 5-6 ms" 1 . 



3.2. Two massive planets in orbit around HD 147018 

HD 147018 has been observed with CORALIE at La Silla 
Obervatory since May 2003. Six radial-velocity measurements 
with a typical signal-to-noise ratio of 25 (per pixel at 550 nm) 
were obtained with CORALIE-98 leading to a mean measurement 
uncertainty of 6.1 ms -1 , including photon noise and calibration 
errors. An additional 105 radial-velocity measurements were ob- 
tained with CORALIE-07 with a mean signal-to-noise ratio of 52 
leading to a mean measurement uncertainty of 3.4 ms _1 . An ex- 
ternal systematic error of 5 ms _I was quadratically added to the 
radial velocity uncertainty before performing the period search 
and the model adjustment. It took more than 1400 days to real- 
ize the importance of this target. Indeed, the first 5 measurements 
only showed a quadratic drift that betrayed the presence of a long 
period companion. We had to wait the sixth measurement, taken 
in June 2007, to realize that we missed a planet due to an inade- 
quate temporal sampling ( compared to the period and the phase 
of the planet). In the following months, the presence of a second 
long period companion, with similar radial velocity amplitude 
was discovered. It took another one and half year to disentangle 
the two orbital solutions and to characterize the second planet 
orbital parameters. 

The planetary system consists of two massive giant planets 
with respective semimajor axis a = 0.24 AU and a - 1.92 AU. 
The first planet is eccentric with e = 0.469 and a period of 
P = 44.24 days. It has a minimum mass m/,. sin (J) = 2.12 Mj up . 
The second planet has a much longer period (P = 1008 days) and 
a minimum mass of m c . sin(/) = 6.56 Mj up . Its orbit is slightly 
eccentric (<? = 0.13) which could betray the presence of interac- 
tions between the two massive planets. 

Figure 2 shows the CORALIE radial velocities and the the ad- 
justed 2 planet-keplerian model. The residuals to the model show 
a level of variation of <x = 7.4 ms _1 , yielding a reduced^ of 1.28. 
The orbital elements for HD 147018 b and HD 147018 c are 
listed in Table 2. Error bars were computed using 5000 Monte 
Carlo simulations and a confidence interval of 68.3%. 

3.3. A long period and massive planet around HD 1 71238 

HD 171238 has been observed since October 2002. Thirty two 
radial-velocity measurements with a typical signal-to-noise ra- 
tio of 17 (per pixel at 550 nm) were obtained with CORALIE-98 
leading to a mean measurement uncertainty of 6.8 ms _1 , includ- 
ing photon noise and calibration errors. An additional 65 radial- 
velocity measurements were obtained with CORALIE-07 with a 
mean signal-to-noise ratio of 47 leading to a mean measurement 
uncertainty of 3.6 ms~'. An instrumental error of 5 ms _1 was 
quadratically added to the radial velocity uncertainty before per- 
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Fig. 2. Radial-velocity measurements of HD 147018 obtained 
with CORALIE-98 (red) and CORALIE-07 (blue). The top panel 
presents the observed radial velocities as a function of Julian 
Date with the best 2-planet-keplerian model (black curve). 
Residuals show a 7.4 ms _1 dispersion. The two bottom figures 
represent the phase folded radial velocities of HD 147018 b (mid- 
dle) and of HD 147018 c (bottom). 



forming the period search and the model adjustment. A clear 
signature is identified in the periodogram at P — 1523 days 
(see Fig. 4) which corresponds to a massive planet (mj. sin(0 = 
2.60 Mj U p) with a semi-major axis of a — 2.54 AU and with an 
eccentric orbit (<? = 0.40). Figure 3 shows the CORALIE radial 
velocities and the corresponding best-fit Keplerian model. The 
orbital elements for HD 171238 b are listed in Table 3. 
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Table 2. Two planet Keplerian orbital solution for HD 147018 
as well as inferred planetary parameters. Confidence intervals 
are computed for a 68.3% confidence level after 5000 Monte 
Carlo iterations. AT is the time interval between the first and last 
measurements, Xr is the reduced^, G.o.F. is the Goodness of Fit 
and cr(0 - C) is the weighted r.m.s. of the residuals around the 
derived solution. The Julian Date is expressed as JD*=JD-2450 
000. C98 stands for CORALIE-98 and C07 for CORALIE-07. 
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The residuals to the adjusted single planet keplerian model 
are however quite large for CORALIE (cr = 10 ms~'), yield- 
ing a reduced x of 1.60. In order to explain such a large dis- 
persion, we have conducted a frequency analysis of the residu- 
als. As can be seen on Fig. 4, a significant amount of energy is 
present around 60 and 160 days. However, no "realistic" keple- 
rian could be adjusted with such periods, discarding the presence 
of additional jovian planets. Furthermore, as explained in sec- 
tion 2, HD 171238 is a DY Drac variable star that could vary on 
time scales of 1 to 120 days. The most likely explanation for the 
large radial velocity dispersion is the presence of stellar spots on 
the surface of the star which is confirmed by a relatively strong 
Call-H re-emission that can be seen in the spectra as illustrated 
in Fig.l 
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Fig. 3. The two diagrams represent the observed radial velocities 
as a function of Julian Date obtained with CORALIE-98 (red) 
and CORALIE-07 (blue) for HD 171238 and HD 204313. The best 
single-planet keplerian model is represented as a black curve and 
residuals show a dispersion of 10 ms _1 for HD 171238 and of 
8 ms- 1 for HD 2043 13. 



3.4. A long period and massive planet around HD 20431 3 

HD 2043 13 has been observed since September 2000. Seventy 
one radial-velocity measurements with a typical signal-to-noise 
ratio of 28 (per pixel at 550 nm) were obtained with CORALIE-98 
leading to a mean measurement uncertainty of 5.2 ms _1 , includ- 
ing photon noise and calibration errors. An additional 26 radial- 
velocity measurements were obtained with CORALIE-07 with a 
mean signal-to-noise ratio of 74 leading to a mean measure- 
ment uncertainty of 3.2 ms" 1 . An instrumental error of 5 ms _1 
was quadratically added to the radial velocity uncertainty before 
performing the period search and the model adjustment. Radial 
velocity measurements betray the presence of a long period 
(P = 1931 days) and massive planet (m&. sin(i) = 4.05 Mj up ) 
orbiting its parent star at a = 3.08 AU. The orbital elements for 
HD 2043 13 b are listed in Table 3. Figure 3 shows the CORALIE 
radial velocities and the corresponding best-fit Keplerian model. 
The residuals to the single planet keplerian model show a level 
of variation (cr = 8.1 ms -1 ), yielding a reduced^ of 1.24. 



4. Concluding discussion 

We have reported in this paper the detection of four extrasolar 
planet candidates discovered with the CORALIE echelle spectro- 
graph mounted on the 1.2-m Euler Swiss telescope at La Silla 
Observatory. 

- HD 147018 b and HD 147018 c are two massive planets part 
of the same system with respective masses mi, sin i = 2.12 Mj up 
and m c sin /=6.56 Mj up . The inner planet has a 44.24 day-period 
and a large eccentricity (e = 0.46) while the outer planet has a 
1008 day -period with a low eccentricity. 

- HD 171238 b is a long period and massive planet with an 
eccentric orbit (P=4.14 years, msini - 2.60 M Jup , e = 0.40). 

- HD 2043 13 b is a long period and massive planet with a low 
eccentricity (P=5.28 years, msini = 4.0 Mj up , e = 0.13). 

It is worth to note that the three parent stars are metal rich, 
which strengthened the case that massive planets tend to form 
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Fig. 4. Periodogram of HD 171238 velocities (a) and of the 
residuals after substraction of the single planet model (b). A sig- 
nificant amount of energy is detected in the residual's close to a 
period of 50 days which is probably activity induced. 

Table 3. Single planet Keplerian orbital solutions for 
HD 171238 and for HD 204313 as well as inferred planetary pa- 
rameters. Confidence intervals are computed for a 68.3% confi- 
dence level after 5000 montecarlo iterations. Ar is the time inter- 
val between the first and last measurements, % r is the reduced 
G.o.F. is the Goodness of Fit and <t(0 — C) is the weighted r.m.s. 
of the residuals around the derived solution. The Julian Date is 
expressed as JD*=JD-2450 000. C98 stands for CORALIE-98 
and C07 for CORALIE-07. 



Parameters 




HD 171238 b 


HD 204313 b 


TC98 


[kms -1 ] 


21.662 ±0.006 


-9.785 ± 0.0018 


7C07 


[kms" 1 ] 


21.641 ±0.002 


-9.762 ±0.0031 


p 


[days] 


1523^0 


1931 ± 18 


K 


[ms- 1 ] 


52.2 ± 1.8 


64.8 ± 1.5 


e 

CO 


[deg] 


400 +00S1 
u.tuo 0065 

47 +9 - 9 


0.131 ±0.023 
-57 ±11 


T 


[JD*] 


5062 ± 20 


3989 ± 62 


a,\ sin; 

/i («) 
m p sin i 
a 


[10 3 AU ] 
[10- 9 M G ] 
[M Jup ] 
[AU] 


6.68 ±0.32 
17.3 ±2.4 
2.60 ±0.15 
2.54 ± 0.06 


11.40 ±0.31 
53.2 ±4.0 
4.05 ±0.17 
3.082 ± 0.055 


N 

Ar 

Xr 

G.o.F 

f (O-C) 


[years] 
[ms- 1 ] 


96 
6.82 
1.60 ±0.07 
7.41 
10.25 


97 
8.06 
1.24 ±0.07 
3.09 
8.11 



around metal rich stars as stated by (Santos et al. (2001, 2005); 
Fischer & Valenti (2005) ). 

On the statistical point of view, the giant planet mass distribu- 
tion decreases with a power law, as illustrated on Fig. 5, with no 
cut-off or change of distribution morphology up to 25 -Jupiter- 
masses. There is no indication in the mass distribution diagram 
(within 6 AU of the parent star) of the presence of the low mass 



tail of a possible brown-dwarf mass distribution. The brown 
dwarf desert is therefore extremely "dry" at the lowest masses 
within 6 AU. However, with only 26 planets candidates with 
masses larger than 7 Mj up , one should be carefull not to over- 
interpret the observations and only a direct measurement of the 
orbital inclination of each planet candidate will reveal their true 
nature. 




m p .sin(i) [Mjup] 

Fig. 5. Frequency distribution of the observed m p .sm(i) val- 
ues for the 117 known exoplanets with masses larger than 2 M Q 
discovered by radial velocities orbiting G&K dwarfs. The un- 
derlying trip, sin (0 distribution, plotted as a red line, is retrieved 
using a non-parametric approach with an Epanechnikov adaptive 
Kernel as described in Jorissen et al. (2001). 

Direct measurement of the orbit inclination of massive planet 
candidates already produced some results. On one hand, Bean 
et al. (2007), with the HST Fine Guidance sensor, conducted a 
set of astrometric measurements of HD 33636 at a fraction of 
milli-arcsecond accuracy. The authors showed that HD 33636 b 
is not in the planet domain and is indeed an M dwarf. On the 
other hand and surprisingly enough, the CoRot space mission 
found a 21.66-Mj up transiting brown dwarf with a 4.26-day pe- 
riod Deleuil et al. (2008). 

However, a systematic monitoring of the massive planet can- 
didates will only be possible with the forthcoming dedicated 
astrometric facilities such as PRIMA (start of operation end of 
2009) (Launhardt et al. 2008) and GAIA (Perryman et al. (2001), 
launch date end of 201 1). Those instruments/telescopes will de- 
termine the real mass of all massive planets candidates as illus- 
trated by figure 6 and provide statistically reliable numbers about 
the distribution of massive planets and of brown dwarfs within 
6 AU. 

Acknowledgements. We thank B. Pernier for participating to CORALIE the 
observations. We are grateful to the Geneva Observatory technical staff, in 
particular to L. Weber, for maintaining the 1.2-m Euler Swiss telescope and 
the CORALIE Echelle spectrograph. We thank the Swiss National Research 
Foundation (FNRS) and the Geneva University for their continuous support 
to our planet search programmes. N.C.S. would like to thank the support 
from Fundacao para a Ciencia e a Tecnologia, Portugal, through programme 
Ciencia 2007. Support from the Fundacao para Ciencia e a Tecnologia (Portugal) 
to P. F. in the form of a scholarship (reference SFRH/BD/2 1502/2005) is grate- 
fully acknowledged. This research has made use of the VizieR catalogue access 
tool operated at CDS, France. 



References 

Bean, J. L., McArthur, B. E., Benedict, G. F, et al. 2007, AJ, 134, 749 
da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609 
Deleuil, M., Deeg, H. J., Alonso, R„ et al. 2008, A&A, 491, 889 



6 



D. Segransan et al.: The CORALIE survey for southern extrasolar planets. 
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